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One-dimensional quantum fluids are conventionally described by using an 
effective hydrodynamic approach known as Luttinger liquid theory. As the 
principal simplification, a generic spectrum of the constituent particles is re- 
placed by a linear one, which leads to a linear hydrodynamic theory. We show 
that to describe the measurable dynamic response functions one needs to take 
into account the nonlinearity of the generic spectrum and thus of the result- 
ing quantum hydrodynamic theory. This nonlinearity leads, for example, to 
a qualitative change in the behavior of the spectral function. The universal 
theory developed in this article is applicable to a wide class of one-dimensional 
fermionic, bosonic, and spin systems. 

The development of the universal effective description of many-body phenomena is a cen- 
tral problem of the condensed matter theory. The hydrodynamic approach known as Luttinger 
liquid (LL) theory C7][2]|5]) is routinely applied to one-dimensional (ID) interacting systems. As 
a crucial simplification, a generic spectrum of the constituent particles is replaced by a linear 
one, leading to a linear hydrodynamic theory, which is nothing but a collection of noninteract- 
ing oscillators. However, to understand a variety of phenomena, such as Coulomb drag between 
quantum wires @, momentum-resolved tunneling of electrons in nanowires ©, and neutron 
scattering off spin chains ©, one needs to take into account the nonlinearity of the spectrum. 
From classical physics, it is known that the existence of nonlinearities may result in qualita- 
tively new phenomena, such as propagation of solitons and appearance of shock waves. These 
phenomena take place in a variety of experimental situations because classical nonlinear hydro- 
dynamics is universal: It is phenomenologically derived from simple assumptions, which do 
not rely on microscopic details. Although description of linear quantum hydrodynamic theory 
requires only quantum mechanical treatment of noninteracting oscillators, formulation of non- 
linear quantum hydrodynamics remains a challenging task because of divergences typical of 
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nonlinear quantum field theories. In this article, we develop a universal theory of ID quantum 
liquids that includes nonlinear hydrodynamic effects, leading to qualitative changes in predic- 
tions for dynamic response functions (e.g. spectral function). 

If the ID quantum many -body problem for fermions is simplified by replacing a generic 
spectrum of particles by a linear one [the Tomonaga-Luttinger (TL) model (00®], it becomes 
solvable at any interaction strength. The Lorentz invariance introduced by this simplification 
protects the existence of well-defined elementary excitations with linear dispersion relation. 
These excitations are quantized waves of density propagating with a velocity v. Adding a 
fermion to a ID system described by the TL model requires creation of multiple elementary 
excitations (|9]IZ9- This can be seen from the form of the fermionic single-particle spectral 
function A(p,uj), which describes the probability of tunneling a fermion with given momen- 
tum p and energy uj into the system [see Supporting Online Material (SOM) for the precise 
definition ([77])]. The spectral function has a power-law singularity at the energy of collective 
excitation uj = vp (see Eq. S4). The corresponding exponent is determined only by the univer- 
sal LL parameter K [the latter is expressed in terms of the density, compressibility, and sound 
velocity v, the three low-energy properties of ID liquid ©]. 

In the phenomenological LL approach, energy scale p 2 /(2m*) is fully dispensed with (the 
effective mass characterizes spectrum nonlinearity at p = and is defined below). The 
conventional justification for such simplification is irrelevance, in the renormalization group 
sense, of the nonlinearity (2). Indeed, the irrelevant terms hardly affect the fermion propagator 
away from the singular lines in space-time, x ± vt \Jt/m„. However, it is the vicinity 
of these lines that defines the nature of singular behavior of the spectral function. We show 
here that for all spinless ID fermionic models with short range interactions the single-particle 
spectral function at p <C kf is universal. In the vicinity of Fermi wave vector +kf and p, uj > 0, 
for example, A(p, uj) is a universal function of a single argument 

A(p,u)<xA(e), e = ^f-^- (1) 

(hereinafter p is measured from the closest Fermi point, and we use units with % = 1). The new 
nontrivial function A(e) is very different from the LL theory predictions, yet it depends only on 
the LL parameter K. The asymptote of A(e) at e ^> 1 does reproduce the LL theory predictions, 
Eq. S4, but at \e ± 1| <C 1 the spectral function is described by power-law asymptotes with 
new exponents. The exponents are different from the predictions of the LL theory but can 
still be analytically expressed in terms of K. We find numerically the universal single-variable 
crossover function A(e), by relating it to the nonlinear dynamics in nonequilibrium Fermi gases 
(T7?lli5||i2ir73]) . We also briefly discuss applications of our results for bosonic and spin systems. 

The spectral function could be measured in tunneling experiments with electrons in nanowires 
© and cold atoms in elongated traps (|76][77|). A closely related object, transverse dynamic spin 
structure factor, is measurable by neutron scattering off ID spin liquids placed in a magnetic 
field ©. The universal crossover function and its analytically obtained asymptotes also pro- 
vide one with a test for numerical methods to evaluate many-body dynamics of ID models, for 
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example, using density-matrix renormalization group algorithms (|7^][i9l) . 

Within a LL approach, fermionic field ^ is expanded by using its components near Fermi 
points as ^(x, t) « ^r(x, t)e lkfX + ^^(x, t)e~ lkfX , and the kinetic energy term in the Hamil- 
tonian is linearized. Solution of the linearized model can be described by using free bosonic 
fields with linear dispersion. Fermionic operators are expressed as exponentials of free bosonic 
fields, and their correlations are easily evaluated. Including the nonlinearity of the spectrum 
of constituent fermions leads to interactions between bosonic fields (|2]|20|). One cannot treat 
such interactions perturbatively in bosonic language in the vicinity of the line cu = vp because 
even in the second order of perturbation theory corrections diverge there (|27|) . Physically this 
happens because conservation laws of energy and momentum are satisfied simultaneously for 
waves with linear dispersion. Thus, two wave packets spend an infinite amount of time near 
each other, leading to an ill-defined perturbation theory. To understand the effects of nonlinear 
spectrum, it is more convenient to work in the fermionic representation. Recently a connection 
between dynamic response functions of ID quantum liquids and well-known Fermi edge singu- 
larity was elucidated (1221I231 ). It allowed one to evaluate dynamic structure factor S(p, cu) and 
spectral function perturbatively in the interaction between fermions. Moreover, it established 
the form of the effective Hamiltonian defining the true low-energy behavior of a liquid com- 
posed of generic particles with nonlinear dispersion relation. For some integrable ID models, it 
is possible to determine the parameters of the effective Hamiltonian nonperturbatively by means 
of Bethe ansatz (1241 1251 Egb . 

The Hamiltonian of the TL model may be recast into the Hamiltonian of free fermionic 
quasiparticles (|27]|2#]|29|) having a linear spectrum: 



Here ^ R ^(x) and ^ R ^(x) are creation and annihilation operators for quasiparticles on the 
right (left) branch, satisfying usual fermionic commutation relations. Colons indicate the nor- 
mal ordering with respect to filled Fermi seas: for right (left) branch all states with negative 



Because the canonical transformation that diagonalizes the TL Hamiltonian is a Bogoliubov 
rotation in the space of particle-hole excitations, such a relation is linear, p R (x) + p L (x) = 
K(p R (x) + p L (x)). Fermionic operators are related to fermionic quasiparticles using "string" 
operators F R(L) (x) as (e.g. for right-movers) 



H 1 = iv / dx [: ^[(x)V^ L (x): -: ^ R (x) W R (x): 



(2) 




*r(z) = Fl{x)¥ R (x), Fl{x) = exp i f dy (5 + p R (y) + 5^p L (y)) 



(3) 



Here we have introduced parameters 



Using Eqs. [3] and H] together with Eq.[2l one can obtain the usual results for Green's function of 
the TL model (|27l) . 

If one wants to consider effects of nonlinearity, one has to include terms that are less rele- 
vant in the renormalization group sense into quasiparticle Hamiltonian. One such term is the 
nonlinearity of the spectrum of quasiparticles: 

H2 = ^-Jdx{ : . (W[)(W L ): +: (W R )(W R ): ) . (5) 

Here m* is the effective mass, which can be related (1201) to low-energy properties as 1/m* = 
v/K l / 2 dv/dh + v 2 / (2K 3 ^ 2 )dK/dh, where h is the chemical potential. 

In principle, there is another term that needs to be included together with Eq. [5j It amounts 
to interaction between quasiparticles created by operators & L R . It can be shown (1281) . however, 
that in the limit of small p interactions between quasiparticles are weak and can be treated per- 
turbatively, along the lines of (122) \23\i . Perturbation theory is valid as long as the interaction 
between the original fermions (created by ^) is short-ranged. Interactions between quasiparti- 
cles are responsible for weak singularities in S(p, uj) near oj = vp ± p 2 / (2m*), large-a; tails of 
S(p, uj) (|22)) . and for possible finite oc p 8 smearing (|23|) of some of the singularities of A(p, uj). 
All these effects vanish as long as one is interested in the scaling limit p — > 0, e — > const, 
see SOM (1771) f° r more detailed discussion. For models with interactions decaying as oc l/x 2 
or slower, non-analytic dependence of interactions on momentum becomes possible, and one 
can not neglect interactions between quasiparticles. This can be already seen from perturbative 
calculations (|23)) . 

The spectral function A(p, uj) gets modified by the spectrum nonlinearity in a profound way 
because the dynamics of the string operators operators F^nJx, t) in Eq.[3]becomes nonlinear. 
Effective mass m* defines the energy scale ~ p 2 / (2m*) near uj = vp where modifications from 
the TL model take place. Because parameters 5± defining F^n) ( x ' 316 universally related 
to K, full form of the crossover written in terms of a variable e is a universal function of K. 
Investigation of the properties of crossover function A(e) is the main subject of the present 
article. 

Before proceeding to discuss the form of the universal crossover, let us consider the main 
new features of A(p, uj) that arise because of nonlinear spectrum. We find that in the vicinity 
of each low-energy region k w (2n + l)kj spectral function A(p, uj) has a power-law behavior 
near frequencies ± [vp ± p 2 / (2m*)] , which is related to orthogonality catastrophe phenomenon 



A(p, uj) oc const + 



OJ 



(6) 



and notations for /i are shown in Fig. [TJ Such power-law behavior results from multiple low- 
energy particle-hole excitations near left and right Fermi points, which are created when "high 
energy" fermion tunnels into the system. 
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•C ■ To understand 



To be specific, let us focus on the vicinity of +kf for p > and to > 0. Because the fermion 
that tunnels into the system has a momentum near +kf and energy of the system increases for 
lj > 0, we need to consider only the correlator (^ R (x, t)\I/ R (0, 0)). 

Let us first discuss the exponent /Z5~+ at the edge u — (vp + 

its origin, one has to understand the states that can be created by \I/ R , when the energy of the 
tunneling fermion is in the vicinity of the edge. From energy and momentum conservation, such 
state is given by a single fermionic quasiparticle with "large" momentum « p and multiple low- 
energy particle-hole excitations with momenta much smaller then p, as indicated in Fig.[2l Then 
one can neglect all other states (|22ll2~3l) and project quasiparticle operators ^ R (x) and ^ L (x) 
onto narrow (of the width much smaller than p) subbands r,d, and I as ^ R (x) ~ Vv(#) + 
e tpx d(x), ^ L (x) 

The effective Hamiltonian determining the evolution of these states is obtained by projecting 
H 1 + H 2 onto subbands r, /, and d and linearizing the corresponding spectra: 

H r j = iv I dx (: ^(x)V^(x): -: jj>l(x)V^ r (x): ) , (7) 



H d = J dxd\x) 



p 2 . / p 



vp H i [ v H V 

2m* V m, 



d(x). (8) 



The Green's function factorizes as oc e ipx (d(x, t)d^(0, 0))jj d {F r (x, t)F}(0, 0)}jv . To obtain 

string operators F r , Fj from Eq.[3l one should keep only r and / components of the density there. 

-■( p 2 \ t 

The free-particle correlator (d(x, t)d^(0, 0)) jj equals oc e V 2m * ) 5[x—(v + ^- )t] , and 
string correlator can be bosonized and evaluated © in a usual way as (F r (x, t)Fj(0, 0)) g r ( \ x= ( v+ je_) t oc 
t -(5_/(27r)) 2 -(5+/(27r)) 2 Ta y n g Fourier transform of (^ R (x, t)^ R (0, 0)), we obtain the universal 
exponent 

— ='-(:- iVfeV- (9) 




Analogously, exponent /i ,+ for lj — [vp — ^— ) <C tt— is determined by configurations 



2m* / 2m* 

with one quasihole with the momentum « — p, two quasiparticles near right Fermi point, and 
low-energy particle-hole excitations. One can again reduce the problem to three- subband model 
and bosonize states near right and left Fermi points. This way, one obtains the exponent 

^-(y 2 -( 2 -y 2< - 3 ' 

New exponents given by Eqs. [9] and [10] are clearly different from the result for the TL model 
in Eq. S4, which corresponds to the exponent 1 — (27r)] 2 . 

Configurations responsible for the remaining exponents 7Jo~, /io,- consist of "high energy" 
particle-hole excitation on the left branch, particle at the right Fermi point, and low-energy 
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excitations on left and right branches. Singularities near k ~ (2n + l)kf also include n low- 
energy particle-hole pairs with momentum rs 2n&/. All exponents can be obtained by using 
projections onto three-subband models, and the results are summarized in Tabled] 

We now discuss the results for the universal crossover function A(e) in the vicinity of +kf 
for p, to > [details of the derivations are available in SOM ([771) ]. The answer is defined by a 
universal function D(y), determined only by 5 + and normalized as j\ D (y) dy = 1. By using 
D (y) , spectral function can be written as a convolution of contributions from the left and right 
branches. Universal function A(e) in Eq.[T]is related to D(y) as 

A{e)= J 1 ^dyD{y)e{e-y){e~y)^ 2 - 1 . (11) 

One can analytically obtain limiting behavior of D(y) for y — > ±1 from Eqs.[9]to[TT]as D(y) oc 
(1 =F y) d± for V ±1; where 

^=(y 2 -^-=( 2 -y 2 - i>3 - < i2) 

At moderate interaction strength, /Z^f" > 0, function A(e) diverges at £ = 1. Then the ratio 
of the prefactors above and below the singular line is universal, 

aq. + |fe|) = r((fe) 2 )r(i-(fe)^) 
w-o^i-i&i) r ((tr)r(i-(t )2 )' 

To evaluate -D(y) away from the edges, one should be able to calculate the dynamics of 
chiral vertex operators (|77l) . For a nonlinear spectrum, this is a very nontrivial problem, the 
analytic solution of which is not known. Similar correlators have attracted attention recently 
(|72~1I7?1) , and their connection to the nonlinear quantum shock wave dynamics and nonlinear 
differential equations has been discussed. Although it might be possible to proceed similarly for 
the evaluation of D(y), it is not clear whether nonlinear differential equations obtained this way 
will have an analytic solution. We use an alternative approach of (|7?]|75|) , which allows us to 
develop a representation of D(y) in terms of certain determinants built of single-particle (rather 
than many -body) states. These determinants can be evaluated numerically, which practically 
solves the problem of finding D(y). Representative results for D(y) and A{e) for K = 4.54 are 
shown in Fig. [3] 

The universal Hamiltonian given by Eqs. [2] and [5] can be also used to describe gapless 
bosonic and spin— | systems away from particle-hole symmetric ground states. We present 
main results on singularities of their dynamic response functions in SOM (|771) . 

We have constructed universal low-energy theory of a wide class of interacting ID quantum 
liquids without resorting to the simplifications of the Tomonaga-Luttinger model accepted in 
the phenomenological Luttinger liquid description. Unlike the latter, we keep the nonlinear 
dispersion relation of the fermions intact. The replacement of the dispersion relation by a linear 
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one, oj = vp, results in an artificial introduction of Lorentz invariance into the system. Although 
not affecting the low-energy behavior of local properties (such as the local tunneling density of 
states), the introduced symmetry alters qualitatively the predictions for the momentum-resolved 
quantities, such as the spectral function. Keeping the nonlinearity allows us to find the generic 
low-energy behavior of the dynamic response functions of a system of interacting fermions, 
bosons, and spins. Possible extensions of our theory should be able to describe the effects 
of finite temperature, spin systems at particle-hole symmetric points, systems with long-range 
interactions, and fermions with spin. 
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Figure 1: Spectral function. (A) Spectral function A(k, uj) in momentum-energy plane. Shaded 
areas indicate the regions where A(k, uj) ^ 0. The region with uj > corresponds to the particle 
part of the spectrum, and the region with uj < corresponds to the hole part of the spectrum. 
(B) Close-up view of the vicinity of k ps +kf, where p = k — kf. Notations of fi indicate 
which exponents presented in Table \T\ should be used in Eq. [6] Notations for exponents near 
k ~ (2n + l)kf are obtained by substituting corresponding n instead of n = 0. 
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Figure 2: Reduction to the effective Hamiltonian. We show excitations contributing to the 
singularity at uo — (vp + ^-J <C The Hamiltonian given by Eqs. [2]and[5]is reduced to 
the three- subband model in Eqs. [7] and [8j 
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Figure 3: Universal crossover. (A) Universal crossover function D(y) for K = 4.54 and the 
corresponding values 5+/(2n) = —0.3 and 5-/(2n) = —0.83; see Eq. HI Exponents d± 
defining the asymptotic behavior at y — > ±1 are given by Eq. [T2l (B) Universal function A(e) 
for K = 4.54. Exponents ~[m^+ and yU ,+ defining the asymptotic behavior at e — > ±1 are given 
by Eqs.[9]and[l0l The ratio of prefactors determining the asymmetry of the singularity at e — 1, 
see Eq. [H equals 2.96 for K = 4.54. 
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Supporting Online Material 

Materials and Methods 1: Definitions of the dynamic correlation functions 

We are interested mostly in the zero-temperature spectral function 

A(k, uo) = --lmG rct (k, uo), (SI) 

7T 

where retarded Green's function G Tet (k, uo) is defined by (SI) 

G rct {k, u) = -ij J dxdte^-^ x (V(x, t)^{0, 0) + &(0, 0)*(x, t))9(t), (S2) 
and dynamic structure factor (DSF) 

S{p,uo) = Jdxdte l(uJt ~ px) (p(x,t)p(0,0)). (S3) 

Here \frf(x, t), ^(x, t) and p(x,t) are fermionic or bosonic creation, annihilation and density 
operators, respectively. Energy uo is measured respective to the chemical potential, so A(k, uo) 
for uo > (uo < 0) describes the response of the system to an addition of an extra particle (hole). 



Materials and Methods 2: Universal crossover 

In this section we present the details of the derivations the universal crossover function A(e) in 
the vicinity of +kf for p, uo > 0. Before proceeding to the case of the nonlinear spectrum, let us 
present the derivation of a conventional result for the Tomonaga-Luttinger model, 

A(p, uo)oc(uo- vp)*( K+ Tc-*)-H(u - vp), (S4) 

which allows for an easy generalization to the nonlinear case. 

Retarded Green's function for fermions near +kf can be written (52) as a product of two 
terms, determined by left and right Fermi points. Due to linear spectrum, they depend on 
combinations v t + x and v t — x, respectively: 

Gf( X ,t) X -i e[M + x) , e{vt ~ X \ . (S5> 

(i(vt + x) + 0)^ 2 {i{vt -x) + 0)(— ) 2+1 



If one defines L(x, t) and R(x, t) as 

1 

(i(vt + x) +0) ( 2^ (i(vt-x) + 0) 



L(x,t)oc -g—,R(x,t)(x -j——, (S6) 
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then imaginary part of Fourier transform of G^(x, t) can be represented as a convolution of 
two Fourier transforms of L(x, t) and R(x, t) : 

ImGg (p, w) = -i f ^L(p -p,u>- u)R(p, u), (S7) 

where real positive functions R(p, u) and L(p, to) are equal to (up to a positive cut-off dependent 
prefactor) 

L(p, a;) oc 5 (a; + vp)9(u> — vp)(u> — vpy , (S8) 
R{j)i uf) oc 5 (uj — vp)9(uo + up) (lo + vpy~' . (S9) 

Using Eq.[ST] spectral function A(p, to) can be written as 

A(p,u) = - I ^L(p-p,u;-Co)R(p,Co). (S10) 

Physically, A(p, to) describes the probability of tunneling of a fermion with total energy u and 
momentum p. Eqs. [S8llS9l then mean that excitations which are created on the right (left) branch 
should lie on the respective mass shell and have positive (negative) momenta. From energy and 
momentum conservation laws, nonzero contribution to A(p, to) for uj > v p comes only from a 
single point in the integral in Eq. lSlOl which correspond to the following arguments of functions 
R(Pk, lur) and L(p L , cul) in the integrand: 

(u + vp) 

^R = vp R = , (Sll) 

~ vp) , c10 v 
UL = -vp L = . (S12) 

Since only L(p L ,to L ) is singular for uj — > vp, for the Tomonaga-Luttinger model only the 
contribution due to the shake-up of low-energy excitations near the left Fermi point controls the 
exponents at to ps vp. For nonlinear spectrum, shake-up contributions from both left and right 
Fermi points determine the exponents, see e.g. Eq. 9. 

For nonlinear spectrum, G R *(x, t) can still be represented as a product of two terms de- 
termined by left and right Fermi points. The primary modification which takes place is that 
delta-functions in Eqs. IS8llS9l get broadened. Indeed, for nonlinear spectrum the total momen- 
tum of several left (right)-moving quasiparticles doesn't completely define their total energy, 
and the latter is allowed to vary up to ±pL(R)/2m*. However, if one is interested in the scaling 
behavior of A(p,u>) for 

p 2 

co-vp = e- — , (S13) 

one can neglect the finite width of function L, since according to Eqs. IS 121 and IS 131 it is of 
the order pi oc p 4 , which vanishes in the scaling limit p — > 0,e — > const. Broadening of 
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function R, on the other hand, is important. As a result of it the momentum and energy on each 
branch can vary on the order ~ p 2 / (2m,*?;) and ~ p 2 / (2m*), respectively, around the values of 
Eqs. lSlllS12l 

To characterize the broadening of delta-function in Eq.[SH let us introduce a dimensionless 
positive function D(y) defined by 



P 



R(p,vp + y—) ccD(y). 



(S14) 



Since for total momentum p on the right branch vp ±p 2 / (2m*) is the highest (lowest) possible 
energy of a set of quasiparticles, D(y) ^ only if y E (—1, 1). Universal function D(y) is 
determined only by 5 + , and we choose it to be normalized as 



f_ i D{y)dy = l. 



(S15) 



We now discuss how to reduce the evaluation of D(y) to a single-particle problem and solve 
it numerically. We use periodic boundary conditions on a circle of length L. Since from now 
on we will be dealing only with fermions at the right branch, we drop index R for clarity of 
notations, and set m* = 1/2. 

The chiral vertex correlation function which determines D(y) can be written as 



R'{x',t) = (e iH2t ^(x')exp 



—i 



dyS+p(y) 



e~ iH2t exp 



i / dy5+p(y) 



* f (0)>, (S16) 



where one has to average over filled Fermi sea on the right branch. In Eq. IS 1 61 we took into 
account the effect of the linear- spectrum Hamiltonian Hi, Eq. 2, by shifting x' = x — vt. 
Universal function D(y) is determined by the Fourier transform of R'(x', t) as 



R'(p,t) = J dx'e- ipx 'R\x',t) oc J 1 ^e ip2yt D{y)dy. 

In momentum space, R'(x', t) can be written as 

R'(x',t) = He^^e^'i^^e-^e^l), 

k<0 p,p' 

where operators A, B and H act in a many-body Hilbert space as 



(S17) 



(SI 8) 



A 



5 + x - 2tt 



p,p 



2ir L(p - p 
2n p ^ p , L{p-p>) p p Z.'WVp" 



p.p 



(S19) 

(S20) 
(S21) 



p,p' 
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while a, b and h are matrices acting in a single-particle Hilbert space. We introduce the density 
matrix 



n 



Z I — Up 



p , (S22) 



where n p is the occupation number of mode p, which we will be set to 9(—p) at the end of the 
calculation. Then Eq. IS18l can be written as a trace over full many-body Hilbert space as 

R>( x ',t) = J] e ik2t Y.e ip ' x 'Tv(^ p ,e B e- mt e A ¥ p p) . 

k<0 p,p' 

We use the relation 

¥ p p = pe x ^ p , (S23) 
which effectively restricts summation to p > 0, and the cyclic property of a trace to get 
R'(x',t) = l[e ik2t £ e ip ' x 'Tr(e B e- mt e A p^ p/ ). 

k<0 p>0,p' 

This trace over many-body Hilbert space can be written via determinants of matrices acting in 
a single-particle Hilbert space as (S3-S7) 

R'(x', t) = J] e ik ' H e ipV (l - np) x 

A:<0 p,p' 
„6 „—iht „a £l\ (~ i „— &Aht „— 6/ 



Det(I -n + e b e- M e a h)(h + e- a e lht e- b (l-n)) , , (S24) 



v ,p 



where n is a diagonal matrix with n p on the diagonal. 

To extract D(y), one needs to evaluate R'(p, t) defined by Eq. lS17l at times 



t a = -f (S25) 



for integer a. According to Eq. lS17[ it corresponds to Fourier series coefficient of D(y) : 

R'(p, t a ) (x D a = e i ^D(y)dy. (S26) 



Function D(y) can be written in terms of D a as 



D(y) = -f + E Re \D a e-^ y \ (S27) 



a=l 



Since for small enough 5+ function D(y) has a singularity given by Eq. 12 for y — > 1, one 
expects 

D a oc (-l) Q «- ( fe )2 , for a -> oo. (S28) 
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We evaluate R'(x', t a ) for various x' using finite-dimensional Hilbert space of the size up 
to ~ 300, and obtain its Fourier transform R'(p, t a ). Due to periodic boundary conditions and 
finite size effects, asymptote given by Eq. lS28l doesn't hold for largest a obtained numerically. 
However, we find an excellent fit for sufficiently large a as 

D a (xe ica (-l) a (sin |) \ (S29) 

Exponent 7 obtained using such fitting procedure equals (§^) 2 with a very high accuracy. In 
Eq. IS291 c C 1 accounts for a possible shift of the frequency, while finite a accounts for 
finite size effects due to finite khj (2n) considered. We remove the latter effects by smoothly 
substituting D a in Eq. lS27l by 

D a oc e iaa (-l) a (^j 7 , (S30) 

for a larger then some intermediate a t <Ca, and keeping numerical results for smaller a. The 
sum in Eq. |S27| with D a given by Eq. |S30| can be written in terms of polylogarithmic functions. 
We sum the contributions to Eq. lS27l coming from large a using polylogarthmic functions, while 
for smaller a we use numerical results. 

The procedure to extract D(y) described above is very robust, and is not sensitive to partic- 
ular choice of parameters at the accuracy of about ~ 2% or 0.02, whichever is larger, for data 
presented in Fig. 3. As an independent check, it reproduces the result D(— 1) ~ very well. 
The correct value of the exponent gL > 3 characterizing the asymptote D(y — > —1) is harder 
to reproduce. 

Materials and Methods 3: Bosonic and spin systems 

The universal Hamiltonian given by Eqs. 2,5 can be also used to describe gapless bosonic and 
spin— I systems away from particle-hole symmetric ground states. The only modification is 
the existence of an additional Jordan- Wigner "string" operator in the expression, in terms of 
fermions, for creation operator of bosons ^ and for spin raising operator S + , respectively. 
Here we will only discuss the singularities. 

For bosons, existence of new singularities in response functions of an integrable Lieb- 
Liniger (S8) model has been pointed out recently (59) and in low-energy regime expressions 
for the exponents in terms of the Luttinger parameter K have been obtained. Same exponents 
can be obtained using the methods of the current article, which demonstrates their universality. 
Exponents in the vicinity of the low energy region k « 2%nkf are summarized in Table ISTl and 
notations of bosonic exponents fi b are indicated in Fig. [ST] 

For spin— I systems, our results apply generally for the following antiferromagnetic (J > 0) 
Hamiltonian in a finite magnetic field h: 

h = jy: s*s? +1 + sfsf +1 + £ v^sfs* -hJ2s?. 

i i>j i 
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Here S* ,y,z are spin— \ operators, and Vi are assumed to decay faster than 1/i 2 , and to be small 
enough so that the system is gapless. We require finite magnetic field, since otherwise due to 
particle-hole symmetry one quite generally has to* = oo. In this case the regime discussed in 
present article disappears, as has been pointed out recently (S10) for an integrable XXZ model. 
While for fermionic systems without a lattice one expects to* > 0, it is not necessarily the 
case for spins on a lattice. This can change the relative position of the singularities compared 
to Fig. 1. For small enough interactions one expects m* > 0(to* < 0) for negative (positive) 
magnetic field h, although for small enough magnetic fields interactions can reverse the sign of 
to*, see e.g. (Sll). 

We will be interested in transverse dynamic spin structure factor, defined by 

S~ + (k,co) =Y^e~ ikj J dte lwt (S-(t)S+ "(0)). (S31) 

It is nonvanishing at low energies in the vicinity of k = tt, as long as the spin chain remains 
gapless (52,572). Generalization of the approach described earlier leads to 



S + (k, uo) oc const + 



1 



UO 



+ ^ 5_ 



-(«l*-l±<£?) 

for | k — it | <C I. Here we have already expressed parameters 5± as functions of K using Eq. 4. 



Materials and Methods 4: Limits of applicability 

In this section we discuss the limitations of and leading corrections to the universal results. 
One regime, when universal results are not applicable has been already pointed out above, and 
corresponds to to* = oo. Such situation generically arises for spin-| system at half-filling, when 
leading correction to spectrum nonlinearity starts from terms oc p 3 . If leading oc p 2 curvature of 
the spectrum is non- vanishing, then our results quite generically apply for 

f- « 1. (S32) 

k f 

We show below that leading corrections to universal results are suppressed in powers of this 
small parameter. To be specific, we consider the modifications of singularities of fermionic 
A(p, uo) for p > 0, uo > in the vicinity of +kf. 

There are two types of terms which modify the universal Hamiltonian. One type of terms 
corresponds to higher order corrections to single-particle spectrum. Such terms merely shift 
the positions of the singularities, but do not change the exponents. Indeed, reduction to three- 
subband model only requires velocity of d— particle to be different from v. Since this happens 
already for leading spectrum nonlinearity oc p 2 , higher order curvature of the spectrum doesn't 
directly affect /J^jT and // ,+- 
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Second type of terms corresponds to irrelevant interactions between fermionic quasiparti- 
cles. One such term, which has the same scaling dimension as spectrum nonlinearity, is given 
by (S11,S13,S14) 

HU = -i~g'J ^(p R [:^V* L :-:V^ L :] - p L [: ^ R W R : -: W R ^ R : ]) , (S33) 

where g' can be related (SI 1, SI 4) to low energy properties similar to 1/m*, and generally these 
quantities are of the same order of magnitude. Effect of such interactions on e.g. JIo~+ can 
be understood using the methods of Refs. (S15,S16). Indeed, after projection to three-subband 
model interactions lead to modification of phase shift 5- of the order 

A<L 9 -^— ~ f- < 1, (S34) 

v d - (-V) k f 

and thus lead to small corrections to /J^jT, /i ,+- Here v d is the velocity of particle d, which 
equals 

v d = v + ^-. (S35) 

m* 

Less relevant interactions between left and right branches lead to even stronger suppressed 
corrections to the exponents. One should note, that presence of finite g' also leads to oc p 8 
smearing of the singularity JT^I for uo > 0, while singularity // ,+, being a singularity at a true 
kinematic border, remains intact (SI 6). 

Another irrelevant interaction term which modifies the exponents in linear order over p/kf 
arises due to momentum dependence of interactions on the same branch, 

H'L = J dpV{p) (pl(p)pl(-p) + Pk( P )pk(-p)) ■ (S36) 

Interaction V(p) should vanish for p — > 0, and have the symmetry property V(p) = V(—p). If 
one assumes that V(p) is regular, then its expansion starts from the term oc p 2 , and correction 
to the phase shift 5+ is of the order 

A5+ ~ ~ < 1. (S37) 

v d - v k f 

However, for interactions that decay as or slower than oc l/x 2 , momentum dependent part of 
V(p) doesn't have to be regular. Indeed, for models with interactions decaying as oc l/x 2 , one 
has V(p) ~ \p\, which leads to finite A5+ in the limit p — > 0, and a finite modification of 
the universal exponents. Thus our universal results do not apply to Haldane-Shastry (S17,S18) 
or Calogero-Sutherland (SI 9) models. In the case of the latter, this can be seen from explicit 
calculations (S16,S20). 
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b 



1 - \ (2„ - 1 - 2„^Y - l (4^)" 
1 _i(2„ + l- 2n fctt) 2 -i(2 

i-K 2 '^ 1 - 2 "^) 2 -^^) 

1 _.( 2 „-l- 2n ^) 2 -l(2 



2 i <5 + -<5-^ 2 



Table SI: Universal exponents for bosonic spectral function. Notations are indicated in Fig. [ST] 
and parameters 5± defined by Eq. 4 are functions of K only. Note that n b n , = //_ n _, which 
follows from the A; — > — A; symmetry. 

Finally, we note that predictions of universal Hamiltonian for S(p, uj) for small p can be 
checked using sum rules. Universal Hamiltonian given by Eqs. 2,5 predicts (527) that S(p, uj) 
at any interaction strength approaches the form characteristic for free fermions 

S(p,uj) = ^d{^- - \u)-v\p\\l (S38) 
p 2m* 

once p becomes small enough. One can check, that for Galilean-invariant systems this result 
explicitly satisfies f-sum rule (522) 

Qf \ duJ 1 ^ 2 n P 2 , C orn 

wS(p,u)— = —vKp = — . (S39) 
2tt 2tt 2m 

In addition, for all systems compressibility sum rule (see e.g. Eq. 7.52 of Ref. (522)) is also 
satisfied: 

hm / -S(p, uj — = — = - — , (S40) 

p^oj uj 2tt 2itv 2 oh 

where in last equation we have used relation of K to compressibility, see e.g. Eq. 2.59 of 
Ref. (52). 
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Figure SI: Bosonic spectral function. (A) Spectral function A(k,u) in momentum-energy 
plane. Shaded areas indicate the regions where A(k, uS) ^ 0. The region with uj > (u < 0) 
corresponds to the particle (hole) part of the spectrum. (B) Close-up view of the vicinity of k ~ 
0, where p = k. Notations of /i indicate which exponents presented in Table ISTl should be used 
in Eq. 6. Notations for exponents near k ~ 2nkf are obtained by substituting corresponding n 
instead of n = 0. 
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